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Abstract 

A new two-component system with cubic nonlinearity: 

m t = bu x + \[m[uv - u x v x )] x - \m(uv x - u x v), 
n t = bv x + \ [n(uv - u x v x )] x + \n(uv x - u x v), 
m = u — u„ 



< 



**XX 7 

n = v- v xx , 



where b is an arbitrary constant, is proposed in this paper. This system is shown integrable 
with its Lax pair, bi-Hamiltonian structure, and infinitely many conservation laws. As a 
special case of b = 0, v = 2, a new Lax pair is given for the standard Camassa-Holm 
equation: m t + m x u + 2mu x — 0,m = u — u xx . In the case of b = 0, the peaked soliton 
(peakon) and multi-peakon solutions are studied. In particular, the two-peakon dynamical 
system is explicitly solved and their collisions are investigated in details. In the case of b ^ 0, 
the weak kink solution is also discussed. 
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1 Introduction 

In recent years, the Camassa-Holm (CH) equation pQ 

mt — bu x + 2mu x + m x u = 0, m = u — u xx , (1) 

where b is an arbitrary constant, derived by Camassa and Holm [T] as a shallow water wave 
model, has attracted much attention in the theory of soliton and integrable system. The CH 
equation was implicitly implied in the work of Fuchssteiner and Fokas on hereditary symmetries 
as a very special case [2j. Since the work of Camassa and Holm [I], more diverse studies on 
this equation have been remarkably developed [3]- [13]. The most interesting feature of the CH 
equation (pQ) is to admit peaked soliton (peakon) solutions in the case of b = 0. A peakon is a 
weak solution in some Sobolev space with corner at its crest. The stability and interaction of 
peakons were discussed in several references |14 j -|18 j . 

In addition to the CH equation, other integrable models with peakon solutions have been 
found [19]-|29]. Among these models, there are two integrable peakon equations with cubic 
nonlinearity, which are 

m t = bu x + [m(u 2 - u 2 x )] x , m = u-u xx , (2) 

and 

mt = u 2 m x + 3uu x m, m = u — u xx . (3) 

Equation ([2]) was proposed independently by Fokas (1995) [Mj, Fuchssteiner (1996) [25], Olver 
and Rosenau (1996) [4J, and Qiao (2006) [26] where it was recognized as a nonlinear water 
wave model. Because of the work by Qiao [26] about its generation from the two-dimensional 
Euler equation, Lax pair, M/W-shaped soliton, and peaked/cusped solitons, equation ([2]) regains 
remarkable attention in the study of peakon and integrable cubic nonlinear system ten years 
after their work [241 [25l d] . Equation ([2]) is the first cubic nonlinear integrable system possessing 
peakon solutions. Recently, the peakon stability of equation ([2]) in the case of b = was worked 
out by Gui, Liu, Olver and Qu [30J. In 2009, Novikov [29] derived another cubic equation, 
which is equation ([2D, from the symmetry approach, and Hone and Wang [28j gave its Lax pair, 
bi-Hamiltonian structure, and peakon solutions. Very recently |31j . we derived the Lax pair, 
bi-Hamiltonian structure, peakons, weak kinks, kink-peakon interactional and smooth soliton 
solutions for the following integrable equation with both quadratic and cubic nonlinearity: 

m t = bu x + ifci [m(u 2 - u 2 x )] x + -k 2 (2mu x + m x u), m = u-u xx , (4) 
where b, k\, and k 2 are three arbitrary constants. 
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It is very interesting for us to study the multi-component integrable generalizations of peakon 
equations. In [U I52~j 133] . the authors proposed the two-component generalizations of the CH 
equation. A two-component integrable extension of the cubic nonlinear equation ([2]) with 6 = 
was introduced in |34| . In this paper, we propose the following two-component system with 
cubic nonlinearity 



rat = bu x + ^[m(uv 



n t = bv x + \ [n(uv - u x v x )] x + \n{uv x - u x v), 
m = u — u xx , 
n = v - v xx , 



(5) 



where 6 is an arbitrary constant. This system is reduced to the CH equation (pQ), the cubic 
CH equation ([2]), and the generalized CH equation as v = 2, v = 2u, and v = k\u + k2, 
respectively. We prove integrability of this system by providing its Lax pair, bi-Hamiltonian 
structure, and infinitely many conservation laws. As a special case of b = 0, v = 2, a new Lax 
pair is provided for the CH equation. In the case of b = 0, we show that this system admits 
the single-peakon of traveling wave solution as well as multi-peakon solutions. In particular, the 
two-peakon dynamic system is explicitly solved and their collisions are investigated in details. In 
the case of b ^ 0, we find that the two-component system ([5]) possesses the weak kink solution. 

The whole paper is organized as follows. In section 2, a Lax pair, bi-Hamiltonian structure 
as well as infinitely many conservation laws of equation ([5|) are presented. As a special case of 
6 = 0, v = 2, a new Lax pair is given for the CH equation. In section 3, the single-peakon, 
multi-peakon, and two-peakon dynamics are discussed for the case of 6 = 0. Section 4 shows 
that equation ([5]) possesses the weak kink solution for the case of 6 ^ 0. Some conclusions and 
open problems are described in section 5. 



2 Lax pair, bi-Hamiltonian structure and conservation laws 



Let us consider a pair of linear spectral problems 





U 



V 




(6) 
(7) 
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where m = u — u xx , n = v — v xx , b is an arbitrary constant, A is a spectral parameter, a = 
yl — A 2 6, and 

A =X~ 2 a + ^(uv - u x v x ) + -(uv x - u x v), 

B = - A _1 (u - au x ) - ^Xm(uv - u x v x ), (8) 

C =X~ 1 (v + av x ) + -Xn(uv — u x v x ). 

The compatibility condition of ([6]) and (J7J) generates 

U t -V x + [U,V} = 0. (9) 

Substituting the expressions of U and V given by © and d7J) into Q, we find that Q is nothing 
but equation ([5]). Hence, ([6]) and ([7]) exactly give the Lax pair of ([5]). 

Remark 1. As we mentioned above, equation ([5]) is reduced to the CH equation ([1]) as 
v = 2. By this reduction, we obtain a new Lax pair of CH equation (pQ) from ([6]) and ([7]): 

1 / A~ 2 a — cm + u x . -A _1 (it-aii x ) + Amu \ 

2 I -2A~ 1 + 2Au -A- 2 Q! + ait - u x )' 



For 6 = (namely, a = yl — X 2 b = 1), that Lax pair reads 

2 ^ 2A 1 J 

if X~ 2 — u + u x — A -1 (u — u x ) + Xmu \ 
2\ v -2A" 1 + 2An -A- 2 + n-^ y' 

which is apparently different from the one given by Camassa and Holm in [I]. In other words, 
here we found a new Lax representation for the Camassa-Holm equation. 



Let 

K=( ° a2 "M 

I i - d 2 o r 

) -1 _ _ - \ (12) 

/ <9to<9 md — md m dmd nd + md 1 n + 2bd \ 
\ dnd~ 1 md + nd~ 1 m + 265 dnd~ 1 nd — nd~ 1 n ) 

By direct calculations, we may verify that J and K are two compatible Hamiltonian operators. 
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Proposition 1 Equation ([5|) can be rewritten as the following bi-Hamiltonian form 

T (5Ht 5Hi\ T (5H 2 5H 2 \ T 

where J and K are given by (f 12$. and 

I r+oo 

Hi = - j (uv + u x v x )dx, 

J-oo (M) 
^ r+oo v ' 

H2 = t / [{u 2 v x + u 2 x v x — 2uu x v)n + 2b(uv x — u x v)]dx. 

J — oo 

Let us now construct conservation laws of equation ([5]). Let ui = where 4>\ and 4>2 are 
determined through equations © and (|7|). From ([6]), one can easily verify that ui satisfies the 
following Riccati equation 



oj x = ——XmuJ 2 + auj — 2^ n - (15) 

Equations © and ([7]) give rise to 

a I 1 
(In = -- + -Amu, (ln^i) t = --(A + flw), (16) 

which yields conservation law of equation ©: 

Pt = F x , (17) 

where 

p = mu, 

1 1 ( 18 ) 

F = X (u — au x )uj — —A 1 (atif — an^fx + nu x — ti^f) + —m(uv — u x v x )oj. 

Usually p and -F are called a conserved density and an associated flux, respectively. Let us 
derive the explicit form of conservation densities for the case of b = 0. To this end, we expand 
uj in terms of negative powers of A as follows: 

oo 

oo = J2uj*~ j - (19) 

3=0 

Substituting (fl9j) into (fl~5j) and equating the coefficients of powers of A, we arrive at 



and the recursion relation for uf 



— 1 



IIIUJO 



LJ-i — U) 



*j,x - 5 m £i+fc=j+i, i,k>i^k , 3 > 1. (21) 
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Inserting f)19[) . ()20|) and f)21 1) into (|18p .we finally get the following infinitely many conserved 
densities and the associated fluxes of equation ([5]) : 

po = yj-mn, F = -V-mn(™ - it^), 

m«j — m x n — 2mn 1 . 1 , N 

pi = - , Fi = --{uv - u x v x + uv x - u x v) + -piiuv - u x v x ), {22) 

2mn 2 2 

pj = mujj, Fj = (u - u x )ujj_ 2 + l}Pj(uv - u x v x ), j > 2, 
where ujj is given by (|20|) and (|2T1) . 



3 Peakon solutions for the case of b = 

Let us suppose that the single peakon solution of ([5]) with b = is of the following form 

u = cie-l x - c *l , v = c 2 e-l x - ci l , (23) 



where two constants c\ and c 2 are to be determined. From the expressions of u and v in (I23p . 
we see that their first order derivatives are discontinuous at x = ct. Thus f)23[) can not be a 
solution of equation ([5]) with b = in the classical sense. However, with the help of distribution 
theory we are able to write out u x , m and v x , n as follows 

u x = —c\sgn(x — ct)e~' x_c *', m = 2c\5(x — ct), 
v x = —C2-sgn{x — ct)e ' s ct \ n = 2c2<5(x — ct). 

Substituting (|23p and (|24p into ([5]) with 6 = and integrating in the distribution sense, one can 
readily see that c% and c 2 should satisfy 

c\C2 = —3c. (25) 

In particular, for c\ = c 2 , we recover the single peakon solution u = ±\/—3ce~\ x ~ ct \ of the cubic 
CH equation © with 6 = [301 13T]. 

Let us now assume the two-peakon solution as follows: 

u = Pl (t) e -\ x - q ^ + p 2 {t)e~\ x - q ^ t) \ v = ri (t)e- lx - qi{t)l +r 2 (t)e- lx - q2{t)l . (26) 

In the sense of distribution, we have 

u x = -pisgn(x - qi)e~^ x ~ qi ^ - p 2 sgn(x - q 2 )e~ ][X ~ q2 \ m = 2piS(x - gi) + 2p 2 S(x - q 2 ), 

-I - I -I - I ( 27 ) 

v x = — r\sgn(x — q±)e |x 911 — r 2 sgn(x — q 2 )e |x 921 , n = 2r\5(x — q\) + 2r 2 5(x — q 2 ). 
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Substituting (|26|) and (|27"|) into © with 6 = and integrating through test functions yield the 
following dynamic system: 

Pi,t = lpi(Pir2 ~ P2n)sgn(q 1 - q 2 )e~\ qi ^ q2 \ 
P2,t = \P2{P2T\ - pxr 2 )sgn(q 2 - qx)e~\ q2 ~ qi \ 

qi,t = 

Q2,t = 



-32 



-\p\r\ - \ (PlT2 +P2Tl) e~ l9l ~ 
-\P2T2 ~ \ {P\T2 +P2Tl) e" l<?2 ~ 911 , 

n,t = -\r\(v\T2 - P2r\)sgn{qi - g 2 )e" l<?1 " 921 , 
r 2 ,t = -\r 2 {p2Ti - Pir 2 )sgn(q2 - qi)e~\ q2 ~ qi \ . 

Guided by the above equations, we may conclude the following relations: 

Pi = Dp 2 , pm = A 1 , p 2 r2 = A 2 , 

where D, A\ and A2 are three arbitrary integration constants. 
If A\ = A 2 , we arrive at the following solution of (|28j) : 



i ^} f5 {D 2 A l -A 1 )s g n(C 1 )e-\ c i\t 
n-i 

P2(t) 

n(t) 
r 2 (t) 
Qi(t) 

92 (*) =qi(t) -Ci, 



pi(t) = Se2D 
Pi 

Pi ' 
Ai 

P2 ' 



^ 1 + 2^ (Z?2yll + Al)e ~ |Cl1 



(28) 



(29) 



(30) 



where B and C\ are two arbitrary non-zero constants. In this case, the collision between two 
peakons will not happen since q2 (t) = qi(t) — C\. 

In particular, as A\ = B = D = 1, C\ = 2, (j30|) is reduced to 



Pi{t)=p 2 {t)=r x {t)=r 2 {t) = \, 



1 



+ e" 



t + 1, 



1 



+ e" 



t - 1. 



Thus, the associated solution of ([5]) with 6 = becomes 

u(x,t) = v(x,t) = e -| a: +(i+ e - 2 )*- 1 l +e -h+(l+ e - 2 ) t + 1 l 



(31) 



which has two peaks, and looks like a M-shape soliton solution [26], but not, and see Figure [T] 
for this M-shape two-peakon solution. 
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As A\ = —B = — D = 1, C\ = 2, the associated solution of ([5]) with 6 = becomes 

u(x, t) = v(x, t) = -e-Mi-^ 2 )*- 1 ! + e -\ x+ ^- e ~ 2 ) t+l \ , (32) 

which has one peak and one trough and looks like N-shape soliton solution. See Figure [2] for 
this N-shape two-peakon (peakon-antipeakon interaction) solution. 

As B = 2D = 1, Ai = Ci = 2, (M becomes 



(33) 



Pi(t) = ~P2(i) = e-i e " 2 *, n(t) = 2ra(t) = 2e§ e ~ 2 ', 
91 (*) = - (I + ^ e " 2 ) * + 92(*) = - Q + ^ e " 2 ) ' - 1 ' 



and the associated solution of Q with 6 = becomes 



u(x,t) = e~l e ~ 2t f e -| <B +(§+i e - a )*-iU2 fi -K(§ + f e_a) ' +J 
v(x,t) = el e ~ 2t f2e-| 1 ' + (l + l e ' 2 )'- 1 l + e -|*+(§+iO*+i 



(34) 



From ([33]) . one can easily see that the amplitudes pi(t) and P2(i) of potential u(x,t) are two 
monotonically decreasing functions of t, while the amplitudes r\(t) and ^(t) of potential v(x,t) 
are two monotonically increasing functions oft. Figures [3] and U] show the profiles of the potentials 
u(x, t) and v(x, t). 
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u(x,t) 



VIX.tl 








-5 






-4 






L 







Figure 3: The two-peakon solution u(x,t) in (f34j) . Figure 4: The two-peakon solution v(x,t) in (f34|) . 
Red line: t = —5; Blue line: t = —2; Brown line: Red line: t = —5; Blue line: t = —2; Brown line: 
t = 0: Green line: £ = 2: Black line: t = 5. t = 0: Green line: t — 2: Black line: t = 5. 



If Ai ^ A 2 , we may obtain the following solution of f|28j) : 
Pl (i) = 5e 2^1-^2) 



Pa(*) 
n(t) 
r 2 (t) 

91 (*) 

92(0 



Pi 
£>' 
Al 
Pi ' 

P2 ' 



(35) 



1 , , 3(A 2 £> a + Ai) 
~3 Alt+ 2D(A 1 -A 2 ) 



3(A 2 D 2 + A 1 ] 
:{ ' 2£>(Ai-A 2 ) 



ssnpi - A 2 )t] ( e -§K^-^)*l _ ^ 



where B is an arbitrary integration constant. Let us study the following special cases of this 
solution. 

Case I. Let A\ = 1, A 2 = 4, B = 1, D = \, then 



Pl (f) = ri(i) = l, P2 =r 2 (t) = 2, 
qi{t) = -\t + 2sgn{t) (e~ - l) , 
[ <8s(t) = -§t + 2ayn(t) (e-l*l - 1) . 



(36) 



The associated two-peakon solution of ([5]) becomes 

l(x,t) = V(x,t) = e -h+|^fln(t)(e-l*l-l)| +2e -h+|t-2 Sg r l ( t )(e-|t|- 1 ) 



(37) 
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u(x,t) 




-15 -10 -5 5 10 15 



Figure 5: The two-peakon solution u(x,t) in (|37l) . 
Red line: t — —5; Blue line: t = —2; Brown line: 
t = (collision); Green line: t = 2; Black line: 4 = 5. 



u(x,t) 





r 
F 






II 


6 























-15 -10 -5 5 10 15 

x 

Figure 6: The two-peakon solution u(x, t) in (|39)) . 
Red line: t — —5; Blue line: t — — 1; Brown line: 
4 = (collision); Green line: t = 1; Black line: t = 5. 



As t < and t is going to 0, the tall peakon with amplitude 2 and at peak position q 2 chases 
after the short peakon with the amplitude 1 and at peak position q\. The two-peakon collides at 
the moment of t = 0. After the collision (t > 0), the peaks separate (the tall peakon surpasses 
the short one) and develop on their own way. See Figure [5] for the detailed development of this 
kind of two-peakon. 

Case 2. Let A\ = 1, A2 = 4, B = 1, D = 1, then we have 



P2(t) 



e* e ", r 2 (t) 



Pi(t) 

n(t) 

gi(t) = -|* + §a^n(t) 
The associated two-peakon solution of ([5]) becomes 

f e -\^+\t-\s g n{t){e-\ t \-l 
\i ( _-\ x+ l t -§ 8 gn(t)(e-W-i)\ 



u(x, t) 
v(x, t) 



e2 



:4el e 

r«— 1*1 



l +e" 
+ 4e" 



(38) 



|s+ff 



+|t-|sgn(t)(e-l' 



1*1-1) 



(39) 



For the potential u(x,t), the two-peakon solution possesses the same amplitude e _ 2 e ; which 
reaches the minimum value at the moment of collision (i = 0). Figure [6] shows the profile of the 
two-peakon dynamics for the potential u(x,t). For the potential v(x,t), the two-peakon solution 
with the amplitudes ea e '*' and 4e2 e 1(1 collides at the moment of t = 0. At this moment, the 
amplitudes attain the maximum value and the two-peakon overlaps into one peakon 5e^e~^, 
which is much higher than other moments. See Figures [7] and [8] for the 2-dimensional and 
3-dimensional graphs of the two-peakon dynamics for the potential v(x,t). 
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v(x.t) 




-10 -5 5 10 



Figure 7: The two-peakon solution v(x, i) in ([39]) . 
Red line: t — —4; Blue line: t = — 1; Brown line: 
t = (collision); Green line: t = 1; Black line: £ = 4. 




Figure 8: 3-dimensional graph for the two-peakon 
solution v(x, t) in ([551) . 



Case 3. Let A 



1, A2 = 4, -B = 1, D = —1, then we have 

' ~P2(t) 



3 e — 1*1 



-1*1 



e 2" , r 2 (t) = — 4e 2 e 
-it-fsgn(t) (e-l'l-1), 



Pl(i) 

n(t) 

Qi(t) 

_ (?2 (t) = -|t-|s( ? n(t)(e 
The associated two-peakon solution of ([5]) becomes 

" u(x,t) = el 6 " 141 L-\*+¥+ls 9 n(t)(e-\n-i)\ _ 
v(x,t) = e-l 6 " 1 ' 1 ( R -\^+¥+^9n(t)(e-W-i 



(40) 



1) 



e -|x+|t+| S sn(t)(e-l t l-l)| 
- 4 e -|^+|*+f^n(i)(e-l t l- 



(41) 



For the potential u(x,t), the peakon-antipeakon collides and vanishes at the moment of t = 0. 
After the collision, the peakon and antipeakon reemerge and separate. For the potential v(x,t), 
the peakon and trough overlap at the moment of t = 0, and then they separate. Figures [9] and 
[10] show the peakon-antipeakon dynamics for the potentials u(x,t) and v(x,t). 
In general, we suppose iV-peakon solution has the following form: 

N N 

u{x,t) = J2Pj(t)e' ]x ~ qj{t)l , v(x,t) =J2 r 3(t) e ~ lx ~ qjm - ( 42 ) 
0=1 j=i 

Substituting ([42]) into ([5]) with 6 = and integrating through test functions, we obtain the 

./V-peakon dynamic system as follows: 

Pj,t = \Pj Effc=iK r fc (sgn(qj - q k ) - sgn( qj - %)) e -l«* - «*H«!i-««l, 
1j,t = \Pj r i ~ 5 £ffc=iPi r fc (1 - S W?? - qi)sgn{qj - q k )) e -l«-*l-l©-»l, (43) 
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Figure 9: Peakon-antipeakon solution u(x, t) in Figure 10: Peakon-antipeakon solution v(x,t) in 
([SI]), Red line: < = -6; Blue line: t = -2; At t = (gTJ). Red line: i = -5; Blue line: t = -2; Brown 
they collide and vanish; Green line: t = 2; Black line: line: i = (collision); Green line: t = 2; Black line: 
< = 6. t = 5. 

It is interesting to study whether the above system is able to be rewritten as an integrable 
Hamiltonian system by introducing a Poisson bracket. We will investigate this in the future. 



4 Weak kink solution for the case of b ^ 

In this section, we show that equation ([5]) with b 7^ possesses a weak kink solution. We assume 
equation ([5]) has the following kink wave solution: 



u 



C lS gn(x - ct) (V |x - Ct| - l) , v = C 2 sgn(x - ct) (V^^l - l) , (44) 



where the constant c is the wave speed, and C\ and C2 are two constants to be determined. In 
fact, if Ci 7^ and C2 7^ 0, the potentials u and v in (|44p are kink wave solutions due to 



lim^+oo u = - lim^^oo u = -C±, 
lim^+oo v = - lim^.oo v = -C 2 - 



(45) 



One may easily check that the first order partial derivatives of (|44p read 

(46) 



u x = -de-\ x - ct \, u t = cC x e-\ x - ct \ 
v x = -C 2 e-\ x ~ ct \, v t = cC 2 e-\ x - ct V 



But, unfortunately, The second and higher order partial derivatives of (|44p do not exist at the 
peak point x = ct. Thus, like the case of peakon solutions, the kink wave solution in the form 
of (|44p should also be understood in the distribution sense, and therefor we call (|44p the weak 
kink solution to equation (0) with b 7^ 0. 
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Figure 11: The weak kink solution at t = 0. Black line: u(x, 0); Blue line: v(x,0). 
Substituting (f4"4"j) and (liUl) into ([5]) and integrating through test functions, we may arrive at 



c = -\b, 
C X C 2 = -b. 



2 (47) 



In the above formula, c = — means that the kink wave speed is exactly —\b. In particular, 
we take b = —2 and C\ = 1, then the corresponding weak kink solution is cast to 

u = sgn(x - t) fe -1 * - * 1 - lV v = 2sgn(x - t) (V 1 * - * 1 - l) . (48) 

See Figure [TT] for the profile of the kink wave solution. 

Usually, multi-peakon solutions take the form of superpositions of the single-peakon solutions. 
But, we need to point out that equation ([5]) with 6^0 does not allow the multi-kink solution 
in the form of the superpositions of single-kink solutions like this: 



N N 

U = 

3=1 3=1 



iv iv 

J2Pj(t)s9n(x - q 3 {t)) ( e -l-^WI - l) , v = Y,r ] {t)sgn{x - q 3 {t)) ( e -l*-*WI - l) .(49) 



In fact, substituting (|49p into © and integrating through test functions, we find that the 
solution assumed in the form (|49p is reduced to zero or single- kink solution (|44p . In paper |31j . 
we proposed the peakon-kink interactional solutions in the form of 

N 

u = P0 (t)sgn(x - q (t)) ( e H*-*>(*>l - l) + ^p^e^-^, 

3=1 

for the cubic CH equation ([2|) (see section 4.2 in [31]). Here, we naturally hope to seek the 
peakon-kink interactional solutions for the two-component cubic equation ([S} with b ^ in the 
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form of 

u = Po (t)sgn(x - q (t)) ( e -l*-*>(*)l - l) + Z$ =1 p j (t)e-\*-*®\, 
v = r (t)sgn(x - q {t)) ( e -l*-*(*)l - l) + J2f=i r ji t ) e ~ lx ~ qj{t)l ■ ° ' 

Unfortunately, by direct calculations, we find that such peakon-kink interactional solutions in 
the form (|50p exist only when the potentials u and v satisfy the reduction condition u = kv, 
which reduces the two-component cubic system © to the cubic CH equation fl2J). In other 
words, we have not yet found multi-kink wave solutions and peakon-kink interactional solutions 
to the two-component cubic CH system (J5J) for the general case: u ^ kv. This is an interesting 
and challenging topic, and we will make a further study elsewhere. 

5 Conclusions and discussions 

In our paper, we present a new two-component integrable system ([5]) with cubic nonlinearity 
and provide its Lax representation, bi-Hamiltonian structure, and infinitely many conservation 
laws. The peakon solutions for the system are derived for the case of b = 0. In particular, 
the two-peakon dynamical system is explicitly solved. The M-shape and N-shape peakons are 
analytically obtained, and the collisions between two peakons are discussed and shown through 
their graphs. We also find that this system admits weak kink solutions for the case of b ^ 0. 
Other topics, such as smooth soliton solutions, cuspons, peakon stability, and algebra-geometric 
solutions, remain to be developed for our new two-component integrable system 
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